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Abstract We discuss the unitary Fermi gas made of dilute and ultracold atoms with an infinite s-wave 
inter-atomic scattering length. First we introduce an efficient Thomas-Fermi-von Weizsackcr density 
functional which describes accurately various static properties of the unitary Fermi gas trapped by an 
external potential. Then, the sound velocity and the collective frequencies of oscillations in a harmonic 
trap are derived from extended superfluid hydrodynamic equations which are the Euler-Lagrange 
equations of a Thomas-Fermi-von Weizsacker action functional. Finally, we show that this amazing 
Fermi gas supports supersonic and subsonic shock waves. 
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1 Introduction 



For interacting fcrmions at very low temperature, far below the Fermi temperature Tp, the effects of 
quantum statistics become very important jl[ . When the densities of the two spin components are 
equal, and when the gas is dilute so that the range of the inter-atomic potential is much smaller than 
the inter-particle distance, then the interaction effects arc described by only one parameter: the s- 
wave scattering length qf 0; 0. The sign of ap determines the character of the gas. Fano-Feshbach 
resonances can be used to change the value and the sign of the scattering length, simply by tuning an 
external magnetic field. At resonance the scattering length ap diverges so that the gas displays a very 
peculiar character, being at the same time dilute and strongly interacting. In this regime all scales 
associated with interactions disappear from the problem and the energy of the system is expected to 
be proportional to that of a non interacting fermions system. This is called the unitary regime 

Recently it has been remarked Q that the unitary Fermi gas at zero temperature can be described 
by the density functional theory. Indeed, different theoretical groups have proposed various density 
functionals. For example, Bulgac and Yu have introduced a superfluid density functional based on 
a Bogoliubov-de Gennes approach to superfluid fermions [3; Papenbrock and Bhattacharyya [|| 
have instead proposed a Kohn-Sham density functional with an effective mass to take into account 
nonlocality effects. Here we adopt the extended Thomas-Fermi (ETF) functional of the unitary Fermi 
gas we have proposed few years ago 0; Q which is a functional of the fermions number density n(r) 
and of its gradient. The total energy in the ETF functional contains a term proportional to the kinetic 
energy of a uniform non interacting gas of fcrmions, plus a gradient correction of the von- Weizsackcr 
form Aa 2 /(8m)(Vn/n) 2 @. This app roach has been adopted for studying the quantum hydrodynamics 
of electrons by March and Tosi (To] , and by Zaremba and Tso UM . In the context of the BCS-BEC 
crossover, the gradient term is quite standard [HI; EH; O E3 Il6t Il7l: [iH : The main advantage of 
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taking such a functional is the fact that, as it depends only on a single function of the coordinates, 
there is no limitation on the number of particles N which it can treat. Other functionals, based on 
single-particle orbitals, require self-consistent calculations with a numerical load increasing with N. 

In this paper we review the last achievements obtained by using our ETF density functional and its 
time-dependent versi on M ; Q . Indeed we have successfully applied this density functional to investigate 
density profiles @; H; I20f7 collective excitations [20[ , Joscphson effect [2l| and shock waves [22[ of the 
unitary Fermi gas. In addition, the collective modes of our EFT density functional have been used 
to study the low-temperature thermodynamics of the unitary Fermi gas (superfluid fraction, first 
sound and second sound) [23j and also the viscosity-entropy ratio of the unitary Fermi gas from zero- 
temperature elementary excitations [24j . 



2 BCS-BEC crossover and the unitarity limit 




In 2002 the BCS-BEC crossover has been observed [25[ with ultracold gases made of fermionic alkali- 
metal atoms. This crossover is obtained by changing with a Fcshbach resonance the s-wave scattering 
length ap of the inter-atomic potential: if ap — > - one gets the BCS regime of weakly- interacting 
Cooper pairs, if uf ±oo there is the unitarity limit of strongly-interacting Cooper pairs, and if 
ap — > + one reaches the BEC regime of bosonic dimers. The many-body Hamiltonian of a two-spin- 
component Fermi system can be written as 

N ± ( & \ 

(1) 

where U (r) is the external confining potential and V(r) is the inter-atomic potential. Here we consider 
AT-f = N±. The inter- atomic potential of a dilute gas can be modelled by a square well potential: 

^)={-„ F °&"^: < 2 > 

where r$ is the effective radius. By varying the depth Vo of the potential one changes the s-wavc 
scattering length 

/ tan(r y/mVo/h)\ 
ap = r 1 r^rr/t — ■ ( 3 ) 



r ^mV /h 

The crossover from a BCS superfluid (ap < 0) to a BEC of molecular pairs (ap > 0) has been 
investigated experimentally Q, and it has been shown that the unitary Fermi gas (|of| = oo) exists 
and is meta-stable. In few words, the unitarity regime of a dilute Fermi gas is characterized by 

r < n~ 1/3 < \a F \ . (4) 

Under these conditions the Fermi gas is called unitary Fermi gas. Ideally, the unitarity limit corresponds 
to 

ro = and ap — ±oo . (5) 

The detection of quantized vortices under rotation [26j has clarified that the unitary Fermi gas is 
superfluid. 

The only length characterizing the uniform unitary Fermi gas is the average distance between 
particles d = ?i~ 1//3 . In this case, from simple dimensional arguments, the ground-state energy per 
volume must be 

V ^52m V ' ^ 5 ' K ' 

with ep Fermi energy of the ideal gas, n = N/V the total density, and £ a universal unknown parameter. 
Monte Carlo calculations and experimental data with dilute and ultracold atoms suggest that 
£ ~ 0.4. 
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3 Extended Thomas-Fermi density functional 



The Thomas-Fermi (TF) energy functional of the unitary Fermi gas in an external potential U(r) 
is given by 



.3 Ti 



£± (37r 2 ) 2/3 n 5/3 (r) + U(r) n(r) 

5 2m 



E = J d 3 r 

with n(r) = n^(r) + nj,(r) total local density. The total number of fermions is 

N = J d 3 r n(r) . 

By minimizing E one finds 

e |l (37r2) 2/3 n2 /3 (r) + c/(r)=/2) 

2m 

with /2 chemical potential of the non uniform system. The TF functional can be extended to cure 
the pathological TF behavior at the surface, namely the fact that the density goes to zero at a finite 
distance from the center of the cloud. We add to the energy density the term 



(7) 



(8) 



(9) 



A 



k 2 {Vnf 
8m n 



A^(Vv^) 2 
2m 



(10) 



Historically, this term was introduced by von Weizsacker Q to treat surface effects in nuclei. Here 
we consider A as a phenomcnological parameter accounting for the increase of kinetic energy due the 
spatial variation of the density. The new energy functional, that is the extended Thomas-Fermi (ETF) 
functional of the unitary Fermi gas, reads 



E = d 3 r 



8m n(r) 5 2m 



By minimizing the ETF energy functional one gets: 



A^V 2 +^(37r 2 ) 2 / 3 n(r) 2/3 
2m 2m 



U(r) 



yn(r) = // y/nfr) 



(11) 



(12) 



This is a sort of stationary 3D nonlinear Schrodingcr equation. In recent papers M ; Ja 120 : l23t |24| we 
have determined the parameters £ and A by fitting recent Monte Carlo results [23 [28|; [29[ for the energy 
of fermions confined in a spherical harmonic trap of frequency u in this regime. The main conclusion 
is that the values 

£ = 0.40 and A = 1/4 (13) 

fit quite well Monte Carlo data of the unitary Fermi gas. Having determined the parameters £ and A, 
one can use our single-orbital density functional to calculate various properties of the trapped unitary 
Fermi gas. Ground-state energies and density profiles have been analyzed in Ref. @; H3|, showing a 
very good agreement with other theoretical aproaches [28|; 0] • 



4 Generalized superfluid hydrodynamics 

Here we analyze the effect of the gradient term on the dynamics of the superfluid unitary Fermi gas. 
At zero temperature the low-energy collective dynamics of this fermionic gas can be described by the 
superfluid equations of inviscid hydrodynamics [2|: 



On 
~dt 



+ V ■ (nv) = , 



(37r^n) 2/3 + U(r) + 



'2m 



= 



(14) 
(15) 



i 



where the velocity v is irrotational: V A v = 0, i.e. 



v(r,i) 



2m 



V0(r,f) , 



(16) 



with 8(r, t) the phase of the condensate wave function of Cooper pairs. It is straightforward to show 
that these equations are the Euler-Lagrange equations of the following TF action functional 
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dtd 3 r 



^0 n+ -L(W) 2 n + ^^(3^ 2 ) 2 / 3 n 5 / 3 + UM ■ 
2 8m 5 2m 



(17) 



Clearly, if the space-time variations of the phase are zero one recover the TF energy functional ([7]). 
From Eqs. (fl"4"]l and (| 15[) one finds the dispersion relation of low-energy collective modes of the 
uniform (U(r) = 0) unitary Fermi gas in the form 



ftcol = ci q , 

where £2 co i is the collective frequency, q is the wave number and 



ci 



VF 



(18) 



(19) 



is the first sound velocity, with vf = y^fpr the Fermi velocity of a noninteracting Fermi gas. 

The simplest extension of Eqs. (fT4"f and (fl"5]l are the equations of extended [HI El irrotational and 
inviscid hydrodynamics: 



^ + V-(nv) 



m ^7 v + V 
at 



h 2 V 2 v^ 



+ ^^-(37T 2 «) 2/3 + I7(r) + ^-V 



0, 
= . 



(20) 
(21) 



These equations include the gradient correction and are the Euler-Lagrange equations of the following 
ETF action functional 



.4 



dt d d r 
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+ (37r 2 ) 2 / 3 n 5 / 3 + U(r)n 
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(22) 



By using Eqs. ((20]) and (|21l) one finds that the dispersion relation of low-energy collective modes of 
the uniform unitary Fermi gas reads 



T + — f hq ^ 2 



£ 2mvF J 



(23) 



where the cubic correction depends on the ratio A/£. 

In the case of spherically-symmetric harmonic confinement 



U(t) 



1 



2 2 

m.uj r 



(24) 



we study numerically the collective modes of the unitary Fermi gas by increasing the number N of 
atoms by means of Eqs. (f2T)]) and ([2"Tj) . As predicted by Y. Castin [3(|, the frequency i?o of the monopolc 
mode (breathing mode) and the frequency fl\ dipole mode (center of mass oscillation) do not depend 
on N: 

J?o = 2lu and i?i = to . (25) 

We find instead that the frequency Q2 of the quadrupole (I = 2) mode depends on N and on the choice 
of the gradient coefficient A. We solve numerically [2(J Eqs. (j2T))) and (f2"Tj) with the initial condition 



n(r,t = 0) 



n os r e 



ie(22 2 -x 2 -K 2 ) 



(26) 
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Fig. 1 Quadrupole frequency J?2 of the unitary Fermi gas with TV atoms under harmonic confinement of 
frequency w. Three different values of the gradient coefficient A. For A = (TF limit): SI2 = y/2ui. Figure 
adapted from Ref. [20| . 



to excite the quadrupole mode, where n gs (r) is the ground-state density profile and e a small parameter, 
The results are shown in Fig. [1] where we plot the quadrupole frequency fli as a function of the 
number N of atoms for three values of the gradient coefficient A. The trend shown in the figure is 
captured by the formula 



a-, 



'2 + 6a 



A 



AT2/3 



1 4. 

^ 2 N 2 / 3 



(27) 



where a = 3(6/5) 3 / 2 (37r 2 ) 2 / 3 £/5. This expression is easily derived by using a time-dependent Gaussian 
variational approach [3l[ in the ETF action functional (|22[) . In the limit N —> 00 it gives the TF result 
Q = v^w @, while in the limit N -t it gives Q = 2uj, which is the quadrupole oscillation frequency 
of non-interacting atoms Q- 



5 Shock waves 

One of the basic problems in physics is how density perturbations propagate through a material. In 
addition to the well-known sound waves, there are shock waves characterized by an abrupt change in 
the density of the medium: they produce, after a transient time, an extremely large density gradient 
(th e shock). Shock waves are ubiquitous and have been studied in many different physical systems 
[321 l33l | . Here we investigate the formation and dynamics of shock waves in the unitary Fermi gas by 
using the zero-temperature equations of generalized supcrfluid hydrodynamics, inspired by the very 
recent observation of nonlinear hydrodynamic waves in the collision between two strongly interacting 
Fermi gas clouds of 6 Li atoms [34| . 

Let us consider the unitary Fermi gas with constant density h with U(r) — 0. Experimentally this 
configuration can be obtained with a very large square- well potential (or a similar external trapping), 
such that in the model one can effectively impose periodic boundary conditions instead of the vanishing 
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Fig. 2 Time evolution of supersonic shock waves. Initial condition with a /If — 18 and r\ — 0.3. The curves 
give the relative density profile p(z) at subsequent frames, where If = y h 2 /(mfp) is the Fermi length and 
lof = (-F /h is the Fermi frequency. Figure adapted from Ref. [2^ |. 



ones. A density variation along the z axis with respect to the uniform configuration n can be experi- 
mentally created by using a bluc-dctuned (bright perturbation) or a red-detuned (dark perturbation) 
laser beam [H; . In practice, we perform the following factorization 



by imposing also 



n(r,t) = nj_(x,y)n\\(z,t) 

n_i(x,y) = n± 
n\\(z,t) = n\\ p(z,t) 



(28) 

(29) 
(30) 



such that 

n(r,t) = np(z,t) (31) 

where h = n_LHi|, and p{z,t) is the relative density, i.e. the localized axial modification with respect 
to the uniform density n. We impose periodic boundary conditions along the z axis, namely p(z = 
L z ,t) — p(z — ~L z ,t), with 2L Z the axial-domain length. We set v(r, t) = (0, 0, v(z, t)) with v(z,t) 
the velocity field such that v(z = L z , t) = v(z = —L z , t). Moreover we impose that the initial localized 
wave packet satisfies the boundary conditions p(z = ±L z ,t = 0) = 1 and v(z = ±L z ,t = 0) = 0. 
Because the dimensional reduction is done assuming the uniformess in x,y directions, we shall consider 
the propagation of a plane wave along the z axis. 

Inserting Eq. ((3"Tj) into Eqs. (j20|) and ([2~Tj) one finds the ID hydrodynamic equations for the axial 
dynamics of the superfluid, given by 



P- 



vp + v p 

Cls(p) 2 



0, 

p' = 0, 



(32) 
(33) 
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where dots denote time derivatives, primes space derivatives, and 

c ls (p) = csp 1 ^ (34) 



is the local sound velocity, with c s = q s (1) = \/(,/3vp the bulk sound velocity, vp = \j2epjm is bulk 

Fermi velocity and ep = ^(^tt 2 "-) 2 ^ 3 the bulk Fermi energy. 

We solve Eqs. (j32|) and ([33)) by using a Crank-Nicolson finite-difference predictor-corrector algo- 
rithm [35l ] with the initial condition given by 

p(z,t = 0) = l + 2r 1 e- z ^ (2 ^ , (35) 

and v(z,t = 0) = 0, where n is the bulk density. In Fig. [2] we plot the time evolution obtained with 

cr /If = 18 and 77 = 0.3, with J Ti 2 / (mep) the Fermi length of the bulk system. The figure displays the 
density profile p(z) at subsequent times. Note the splitting on the initial bright wave packet into two 
bright travelling waves moving in opposite directions. There is a deformation of the two waves with 
the formation of a quasi-horizontal shock-wave front. Eventually, this front spreads into wave ripples. 
We have carefully checked that these ripples arc not an artefact of the numerical scheme. Notice that 
before the shock both amplitude and velocity of the two maxima of the two waves are practically 
constant during time evolution. In particular, as we have recently shown [22| that the amplitude of the 
extrema can be written A(rj) = 1 + 77 while the velocity of the extrema reads 

V(r ] ) = ±c s (A{\ + r 1 ) 1 ^ -3) . (36) 

Taking 77 = the velocity of the impulse extrema reduces to the sound velocity: V(0) = c s = yj^/"ivp. 
Moreover, bright perturbations (77 > 0) move faster than dark ones (77 < 0), and the Mach number 
M = V(j])/V(0) of these perturbations in the unitary Fermi gas is simply 

j\/ = 4(l + 7?) 1/3 -3. (37) 

For M > 1, which means 77 > (bright perturbation), one has supersonic waves, while for < M < 1, 
which means i] < (dark perturbation), one has subsonic waves. In the upper panel of Fig. [3] we plot 
the Mach number M as function of the amplitude 77 of the perturbation. Note that since 2rj is the 
amplitude of the initial condition, see Eq. (|35[) . the region 77 < —1/2 is unphysical. 

Let us consider a bright perturbation (77 > 0) moving to the right. The speed of impulse maximum 
V(rj) is bigger than the speed of its tails V(0). As a result the impulse self-steepens in the direction of 
propagation and a shock wave front takes place. The breaking-time T s required for such a process can 
be estimated as follows: the shock wave front appears when the distance difference traveled by lower 
and upper impulse parts is equal to the impulse half-width a/2. This criterion gives [22| 

Ts = 2c s ((l+^/3-l) ' (38) 

In the lower panel of Fig. [3] we plot the period T s as a function of the amplitude 77 of the perturbation. 
The figure shows that as 77 goes to zero the period T s goes to infinity; in fact, in this limit the shock 
wave reduces to a sonic wave (sound wave) which does not produce a shock. 



6 Conclusions 

We have shown that the ETF functional of the unitary Fermi gas can be used to determine ground-state 
properties of the system in a generic external potential U(r). Moreover, the time-dependent version 
of this EFT function, that is the generalized hydrodynamics equations, can be applied to calculate 
sound waves, collective modes and also shock waves. Our generalized superfluid hydrodynamics is re- 
liable when the effects of the temperature T are negligible, namely if T -C T c , with T c the critical 
temperature of the supcrfluid-normal phase transition (T c ~ 0.2 Tp, and Tp — 10 -7 Kelvin for dilute 
alkali- metal atoms). As previously discussed, recently an observation of nonlinear hydrodynamic waves 
has been reported in collisions between two strongly interacting Fermi gas clouds of 6 Li atoms [3~Ij . 
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Fig. 3 Upper panel: Mach number M — v max /c s as a function of the amplitude 77 of the perturbation (solid 
line). Lower panel: period T s of formation (breaking time) of the shock- wave front as a function of the amplitude 

77 of the perturbation. T s is in units of u/c 3 , where a is the width of the perturbation and c s = W^/Svf is the 
bulk speed of sound, with vf the Fermi velocity. Figure adapted from Ref. [22| . 



The experiment shows the formation of density gradients, which are nicely reproduced by hydrody- 
namic equations with a phenomenological viscous term [3J]. Nevertheless, the role of dissipation is 
questionable in this case since the ultracold unitary Fermi gas is known to be an example of an almost 
perfect fluid (36| . We plan to simulate the collision between two strongly interacting Fermi-gas clouds 
of 6 Li atoms by using our single-orbital time-dependent density functional. In particular, we aim to 
demonstrate that it is possible to reproduce the recent experimental results without the inclusion of 
a phenomenological dissipativc term. Another interesting related problem is the determination of the 
surface tension of the unitary Fermi superfluid. We shall consider a semi-infinite domain, derive the 
density profile of the system by using our density functional, and obtain the surface tension as the 
grand potential energy difference between the actual configuration and the uniform asymptotic one. 
Finally, we shall compare our result with previous determinations based on microscopic calculations 



of the normal-supcrfluid interface in population- imbalanccd Fermi gases 37 1 
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